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The Bloch-Horowitz (BH) equation has been successfully applied to calculating the binding 



. energies of the deuteron and ^H/^He systems. For the three-body systems, BH was found to be 

perturbative for certain choices of the harmonic oscillator (HO) parameter b. We extend upon this 
Q ' work by applying this formalism to the alpha particle and certain five-, six-, and seven-body nuclei 

C ■ in the p-shell. Furthermore, we use only the leading order BH term and work in the smallest allowed 



included-spaces for each few-body system (Ohuj and 2huj). We show how to calculate A-body matrix 



H ' elements within this formalism. Stationary solutions are found for all nuclei investigated within 

this work. The calculated binding energy of the alpha particle differs by ~ 1 MeV from Faddeev- 
Yakubovsky calculations. However, calculated energies of p-shell nuclei are underbound, leaving 
p-shell nuclei that are susceptible to cluster breakup. Furthermore, convergence is suspect when 
the size of the included-space is increased. We attribute this undesirable behavior to lack of a 
sufficiently binding mean-field. 
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I. INTRODUCTION 



Difficulties intrinsic to the many-body problem stem from the exponentially-increasing 
complexity of the Hilbert space as the size of the few-body system increases. This in- 
creasing complexity is very daunting yet much progress has been made in developing 
methods that illuminate the nuclear structure of few-body systems. A particular method 
that warrants mentioning is the Green's function Monte Carlo (GFMC) method (see, for 
instance, Refs. 0,0,0] and references within). Using bare 'realistic' phenomenological two- 
and three-body forces, GFMC calculations for certain p-shell nuclei have produced nuclear 
binding energies to a precision better than 2%. This precision, coupled with the GFMC's 
success in reproducing observed nuclear few-body data, has placed the GFMC method as 
the benchmark which other methods use as comparison. Yet even with the consistent in- 
creases in computer resources and speed, GFMC calculations are very computationally tax- 
ing. GFMC calculations on sd-shell nuclei, if possible, are still many years away. Other 
'infinite Hilbert-space' methods succumb to the same problems. 

These difficulties have motivated the use of methods that employ effective interactions 
within severely truncated Hilbert spaces, or included-spaces, thereby alleviating many of the 
complexities of the many-body system. A classic example is the traditional nuclear shell 
model (TNSM). Unfortunately, calculating the effective interaction exactly is as difficult a 
problem as the original 'infinite Hilbert-space' problem. Hence some form of approximation 
must be done to the effective interaction. For example, in TNSM calculations the effec- 
tive interaction is approximated only at the two-body level (the exact effective interaction 
is an A-body interaction, where A is the number of nucleons). Its matrix elements are 
typically ffited to existing nuclear data using a G-matrix0, 0] as a starting point. An alter- 
native exarnple (aiming at a convergence to the exact solution) is the No-Core shell model 
(NCSM) 0, 0, 1^ . In NCSM ab initio calculations, the effective interaction is derived from the 
bare nucleon-nucleon (NN) interaction by performing a unitary decoupling transformation. 
NCSM calculations approximate the exact A-body effective interaction by retaining only the 
two- or two- plus three-body effective terms. Both methods have had resounding success in 
describing nuclear few-body systems. TNSM by and large remains the most popular method 
for working with nuclei in the sd-shell, whereas NCSM results for p-shell nuclei agree very 
well with GFMC. 
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In a recent paper by two of us, we applied an alternative effective interaction given by the 
Bloch- Horowitz (BH) equation[2] to the nuclear three-body svstemjl^. We found that by 
performing a certain ordering of diagrams, the BH equation could be applied perturbatively 
within a finite range of the harmonic oscillator (HO) parameter b. In this paper, we extend 
upon this work by applying the same formalism to the alpha particle and certain p-shell 
nuclei. Unfortunately, the next-to-leading-order and higher-order terms in the BH expan- 
sion are increasingly difficult in p-shell nuclei due to the various many-particle permutation 
operators that come into play. Hence we apply only the leading order BH term to these 
nuclei. Also, we work within the smallest allowed included-spaces so as to alleviate as much 
as possible the complexity of the few-body Hilbert spaces. 

As mentioned in Ref. Il0|, the leading order BH term resembles in some ways that of 
the G-matrix effective interaction!^, |5|]. Yet there are subtleties that differentiate the two 
expressions. In this paper, we further elaborate upon these differences, as well as make 
comparisons to the 'multi- valued' G-matrix of Ref. We show that the leading order 

BH term is a natural extension of the G-matrix in which spectator dependence and cluster 
recoil is explicitly taken into account. These new physical inputs are the main motivation 
for using the leading-order BH term over the G-matrix. Yet their inclusion presents its own 
difficulties. In particular, the leading order BH expression is an A-body interaction. 

In Sec. |H] we give a cursory review of the Bloch-Horowitz equation. We describe our 
methodology and motivate our approximation of the BH equation to its leading-order ex- 
pression in Sec. IIIII In doing so, we compare the BH leading term to the traditional G-matrix 
as well as the 'multi- valued' G-matrix. Since our BH term represents an A-body interaction. 
Sec. IIVI describes how to calculate A-body matrix elements using this term. In particular, we 
show that the degrees of freedom over spectator nucleons can be integrated out completely, 
leaving an average residual two-body interaction. Section |3 presents our results for the 
alpha particle (as well as the three-body system) and certain p-shell nuclei. There are issues 
with regards to convergence of our solutions when we increase our included-space sizes. We 
discuss these issues in Sec. I VII and compare differences between our calculations and NCSM 
calculations. We give concluding remarks in Sec. IVIII To keep the presentation concise, we 
reserve formal derivations for the appendices. 
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II. BLOCH-HOROWITZ EQUATION 



The Bloch-Horowitz equation is 

HeffiE) = P(^H + H^^QH^ P, (1) 

where 

H = t^m, + V,,) = E + - (2) 

i<j i i<j 

Here, the intrinsic Hamihonian H is obtained by summing the relative kinetic and potential 
energy operators Tij and Vij, respectively, over all nucleon pairs i and j. The relative kinetic 
energy is found by subtracting the center-of-mass (cm.) kinetic energy from the total single 
particle kinetic energies. Ma is the total mass of the A-body system. 

In Eq. ^ the operators P and Q = 1—P are projection operators that define the included- 
and excluded-spaces, respectively. They are intimately related to the choice of basis states 
in which calculations are performed. For this work we chose the Jacobi HO basis. This basis, 
defined solely within the relative coordinates, avoids the overcompleteness characterized by 
the the cm. motion intrinsic to traditional slater-determinant HO bases. By avoiding this 
cm. motion, our original A-body problem becomes an (A-l)-body problem. The Jacobi 
included-space (P) in which we perform our calculations is thus defined in terms of all 
relative-coordinate configurations with HO energy quanta Nhu < Ahu. Hence A represents 
the demarcation between included- and excluded-spaces. 

Though many of the complexities of the few-body problem are alleviated by working in 
small included-spaces in which the BH equation is defined, calculating matrix elements of 
Heff is still a very daunting task, if not impossible. The energy E appearing in Eq. ^ for 
example, represents the exact eigenvalue of the system. Obviously, this eigenvalue is not 
known a priori, hence the BH equation must be solved self-consistently. Furthermore, H^/f 
represents an A-body interaction, regardless of the fact that the original bare H may have 
been an (n<A)-body interaction. The most difficult aspect of the BH equation deals with the 
resolvent appearing in Eq. ^ [1/{E — QH)]Q. This resolvent represents the full many-body 
interacting Green's function. The Q operators appearing in this expression implies that 
this propagator resides in the excluded-space, which is of infinite size. Clearly solving for 
this operator is as difficult as the original 'infinite-Hilbert space' few-body problem. Hence 
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some form of approximation of the BH equation is in order. The next section shows how we 
appproximate the BH equation to a more tractable ladder summation. 



III. APPROXIMATING THE BLOCH-HOROWITZ EQUATION 

In Refs. jlO. 12 1 we presented arguments for using a particular rearrangement of Eq. ^ 

where 

= T + T^QT, (4) 
= + 

We stress that this particular form of the effective interaction comes from a formal manipula- 
tion of the BH equation. At this stage no approximations have been introduced. Equation El 
is convenient to work with since it separates out an effective kinetic term and an effective 
potential term. As argued in Refs.[lo[ll^. the operator [E/{E — QT)]P that sandwiches T^^ff 
and Veff can be used to define new basis states that have the correct exponential asymp- 
totic behavior, as opposed to the Gaussian decay of HO wavef unctions. This corrects for 
HO 'overbinding'. Furthermore, matrix elements of T^ff are relatively easy to calculate. In 
Ref. [lol| we gave the analytic expression for calculating three-body matrix elements of T^ff. 
In Appendix El we generalize this expression to A-body systems. 

The crux of the problem deals with K//, since the same troublesome resolvent described 
in the previous section appears in this operator. For the three-body system of Ref.Q], we 
invoked Faddeev decompositions to generate a perturbative expansion of the K//. In what 
follows, we show an alternative approximation to V^ff which is completely equivalent to the 
leading-order term of the expansion of Ref.^^, but generalized to A-body systems. 

A. Retaining the ladder summation 

We begin by expanding the resolvent in powers of QV, 
^ Q = ^ ^ 0+ ^ ^^^ QV ^ ^ 0+ ^ ^^^ QV ^ ^^^ QV ^ ^ + ... 



E-QH^ E-QT^ E-QT^ E - QT^ E - QT^ E - QT^ E - QT 

(6) 



We retain from this expansion only the terms that contribute to a ladder summation, as 
shown diagrammatically in Fig. ^ Here only particles 1 and 2 are interacting via repeated in- 
sertions of the potential V12, while particles 3 through A are spectator nucleons. This infinite 
ladder sum is collectively grouped into the expression Vi2,e//- Hence Veff is approximated 

as 

A 

Veff - ^^='^ff (7) 



A{A - 1) 



^12,e//, (8) 



2 

where the A{A — l)/2 factor comes from using anti-symmetrised states. Obviously the choice 
of ij = 12 in Eq. |H1 is arbitrary. 

There are three distinct ways in which to do the ladder summation. The first method 
assumes that the spectator particles are completely irrelevant. They do not contribute in 
any way to the ladder interactions between particles 1 and 2. Such an approximation is akin 
to the G-matrix interaction of traditional shell-model calculations. The resulting interaction 
is completely two-body in nature which, from a computational standpoint, greatly simplifies 
the problem. Furthermore, the effective interaction is much softer than the original bare A^A^ 
interaction, thereby reducing the coupling of the high momentum matrix elements to the 
included space. This softening of the core is a generic property of the ladder sum. Indeed, 
it can be shown that in the limit of an infinite hard core, matrix elements of this ladder 

Vl2,eff = ^ rj. + ^ rp Q12V12— ^ Ql2 + 

^ rj, <5l2^12^ ^ rj. <5l2Vl2— ^ rj. Ql2 + • • • (9) 

— (^12il2 -£^ — (^12^12 -C/ — VI2J12 

Note the replacement of Q with Q12 in the previous equation. Since the presence of spectator 
nucleons is completely ignored, the many-body projection operator Q acts only within the 



We label 



two-body sector. This condition is valid only in the A = excitation spaces 
this operator as Qi2- Since this interaction is only two-body, it is a poor approximation to 
the original -ffe//, since the latter is an A-body interaction. Furthermore, since T is replaced 
by T12, the energy E can no longer be associated with the energy of the A-body system. 
Historically, the missing many-body correlations have been 'mocked' up by introducing 
fictitious parameters, such as starting energies uj. These further approximations eliminate 



any self-consistency, as well as manifest A-body anti-symmetry. Attempts at improving 
the effective interaction via perturbation have had mixed success as well. Barrett and 
Kirson[^ have shown that matrix elements of third order (in G) diagrams can be as large 
as second order diagrams. These inconsistencies have led shell model practitioners on a 
more phenomenological track, using effective two-body matrix elements derived by fitting 
to nuclear spectra. Clearly any underlying connection to the original bare A^A^ interaction 
has been lost at this point. 

Since the included space is defined in terms of many-body states (and not just two- 
body states) with total HO quanta below A, the actual HO configuration of the spectator 
nucleons restricts the possible configuration in which the two-nucleon cluster may interact. 
This is especially true for multi- valued excitation spaces {i.e. A > 0). For example, within a 
Ahuj = ATioj included space, a 2huj spectator state restricts the possible two-body interacting 
cluster to either O^a; or 2huj, since the sum of the HO quanta of the spectator and interacting 
clusters must add to a value < A = 4. The multiplicity of the two-nucleon state implies that 
the G-matrix is multivalued. This implicit spectator dependence can be included within the 
ladder sum by replacing Q12 by the full many-body Q within Eq. IHl 

^^QVr^^r^Qy^^^^Q^- (10) 

Equation represents the second method for performing the ladder summation. Zheng 
et a/.[ll| have shown that there are significant improvements to using the multi- valued G- 
matrix over its single-valued cousin. In particular, the calculated low energy excitation 
spectrum of certain few-body systems were lowered relative to the ground state when using 
the multi- valued G, bringing energies into better agreement with experiment. Despite its 
spectator dependence, however, Eq.^J still represents a two-body operator since it is diag- 
onal in all spectator quantum numbers. Hence induced multi-particle correlations are still 
missing. 

In the center of mass frame, the relative kinetic energies of the interacting two-nucleon 
cluster and spectator cluster are correlated. This is especially true for light systems. This 
cluster recoil can be incorporated into the ladder sum by using the full kinetic energy operator 
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T in Eq. HOI 

Note that since the full Q and T are present, Eq. ^2 represents a multi-valued, A-body 
operator. The A-body nature of this operator implies that some induced multi-particle cor- 
relations are being incorporated at this level of approximation. However, these correlations 
do not come from potential interactions; rather, they stem solely from the fact that the 
kinetic energies of the clusters are correlated. This is the third and final approximation to 
Veff (Eq. 13 which can be obtained from a ladder summation. It is also the approximation 
of Veff that we investigate in this paper. 



IV. CALCULATING A-BODY MATRIX ELEMENTS OF Fi2,e// 



The infinite ladder sum in Eq. ^2 can be expressed in closed form 



Vl2,eff — tl2 — tuGoP— — - — PGoti2, (12) 

i + i oo 



where 



Go = (13) 
To = PGoP, (14) 

tl2 = Vi2 + Vi2Goh2, (15) 

Too = PGotuGoP. (16) 

This is precisely the leading-order BH term of Ref.Q]. Here Go represents the free many- 
body propagator, while Tq consists of included-space overlaps of Gq. The t-matrix for free 
particles satisfies the Lippman-Schwinger equation that depends on V12 and the many-body 
propagator Gq, and Too represents included-space overlaps of ti2 weighted by two insertions 
of Go. 

The many-body nature of the operator Vu^eff comes directly from the included-space 
overlaps of Goti2Go. In momentum space ti2 is not only a function of the interacting two- 
nucleon cluster's initial and final momenta, p and p' respectively, but also of the individual 
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spectator Jacobi momenta qi 



A-2 



tu{p,p;E-J2i/'^Mn). 



(17) 



i=l 



Here M„ is the nucleon mass. This dependence on spectator momenta imphes that overlaps 
of this operator are not diagonal in all the spectator's quantum numbers, and in particular, 
their principle quantum numbers can change by arbitrary amounts. In principle, calculat- 
ing included-space many-body matrix elements of this operator would require numerically 
evaluating A-dimensional nested integrals in momentum space (the same is true in coordi- 
nate space). Such integrals can be done via Monte Carlo methods, though these methods 
would be very computationally taxing and probably inaccurate due to the oscillatory na- 
ture of the HO wavefunctions. In Appendix ^ we show an alternative method in which 
the multi-dimensional integrals over the spectator nucleons can be reduced to a sum over 
one-dimensional integrals. The result is valid regardless of the number of spectator nucleons, 
and is given by 

J dqi--- dqA-2 qf ■ ■ ■ q\-2 Rn[h (^O^mh (?l) " " " Rn'^_jA-2 {QA~2)RnA-2lA-2 iQA-2) 

A-2 2 

i=l 

C{mi,m[, ■ ■ ■ , mA-2, 



X 



where b is the oscillator parameter and 



;i8) 



C{mi,m[, ■ ■ ■ , mA-2, m'A-2) 



'A-2 



n 



X (_i)S>K+K) 



^2T{n, + l)r(n, + k + 3/2)) 
\ V{ni-mi + l)V{li + mi + ?>/2)V{m, + l) 

v/2r(n^ + l)r(n^ + /, + 3/2")) 
\ r«-m^ + l)r(/i-fm^ + 3/2)r(m^ + l) 



n 



X 



A-2 

n 



r(/, + mi + m[ + 3/2)r(^;.i;(/,- + m,- + m'^ + + 



(19) 



.^2 2r(/, + mi + m[ + E]Jiilj + rrij + m'^ + 1) + ^) 
In Eq. ^Jthe first two products run over all spectator nucleons {i.e. i=l to A-2). Equa- 
tion [THl utilized the fact that matrix elements of ti2 are diagonal in the spectator's orbital 
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quantum numbers k (see Appendix |B}, but not their principle quantum numbers. The 
right hand side (RHS) of Eq. ITHl represents a significant simphfication to the original multi- 
dimensional integrals, since the one- dimensional integrals of the RHS can be done quickly 
and accurately. The coefficients C{mi, m[, - ■ ■ , mA-2, ^'a-t) known analytically and can 
be pre-computed and stored. 

A. Dependence of ti2 on spectator nucleons 

Equation is only a function of the initial and final momenta of the interacting two- 
nucleon cluster [i.e. nucleons one and two). However, its functional form depends implicitly 
on the number of spectator nucleons and their initial HO momentum distribution. In some 
sense, ti2 represents the average interaction of particles one and two under the influence of 
the spectator's momentum distributions. It is not a mean-field interaction, since ti2 is not 
derived from any interactions with the spectator nucleons. We will discuss this point in 
more detail later. 

Figure |21 shows this spectator dependence of the operator 
bhujp'^Golp'; E)ti2{p' ,p] E)Gq{p; E)p'^ , which is dimensionless. Note that uj is the HO 
frequency, defined by fnu = /M^b'^, where Mn is the nucleon mass. We plot this quantity 
rather than ti2{p',p;E) since it is more relevant for our calculations and more visually 
instructive. All plots within Fig. El are for the case when tu is in the ^5*0 channel, plotted 
against initial and final dimensionless momenta bp and bp', respectively. Furthermore, all 
plots were calculated with E=-20 MeV and 6=1.2 fm. Plot (a) shows this operator when 
there are no spectator nucleons present. In this case tu is the standard two-body t-matrix. 
Plot (b) shows how this operator is modified after integrating out the presence of two 
spectator nucleon in the lowest allowed relative Jacobi configuration, the Os state. Plot 
(c) shows the effect of integrating out three spectator nucleons, two of which are in Os 
states and one in a Op state. Finally, plot (d) shows how the interaction is altered with the 
presence if four spectator nucleons are integrated out. In this case, two of the nucleons are 
in relative Os states, and two are in the relative Op states. 

The physical interpretations of the troughs and peaks seen in Fig.|2lare difficult to explain 
due to the complicated structure of the operator Go{p'] E)ti2{p', p; E)Gq{p; E). What is clear 
is that the inclusion of additional spectator nucleons softens this operator. This is consistent 
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with the fact that the possible configurations of the interacting two-nucleon cluster are 
reduced with increasing numbers of spectator nucleons. 



V. RESULTS 

In all our calculations we use the Argonne vis and v'g potentials'^ (Avl8 and Av8') for 
our A^A^ interactions. Furthermore, we ignore isospin symmetry breaking effects included 
in the Avl8 potential since this is a small effect. Table IJ shows the dimensions of the 
A = 0, 2 included-spaces for the ground state quantum numbers of different few-body systems 
considered in this paper. Since the dimensions are very small, our calculations do not demand 
extravagant computing resources and can be done very quickly. To construct our included- 
space anti-symmetrised HO Jacobi wavef unctions, we use a modified version of the routine 
manyeff outlined in Ref.j^]. We make our comparisons with GFMC and NCSM calculations, 
and show experimental results solely as reference points. 



A. S-shell nuclei 

Figure El shows our results of using Vu^eff on s-shell nuclei. All calculations were done 
in A = included spaces and with the AvlS potential. Our results are displayed by points, 
while the 'exact' theoretical results are taken from Ref. [igIi and displayed as solid lines. For 
the three-body system, the upper and lower set of points refer to ^He and ^H, respectively. 
Note that for ^He, the plotted range of b values is much smaller. 

The system is special since in this two-body system there are no spectator nucleons. 
The ladder sum of Fig. ^ is exact in this case, meaning that ^i2,e// is not an approximation 
but corresponds to the exact effective interaction. Thus any calculations of the two-body 
system cannot depend on b, as this is an intrinsic parameter of the HO basis and has nothing 
to do with the A^A^ interaction. Indeed, all our ^H calculations lie on the exact theoretical 
e„e.,y0, regarCe. of the oscUl.c. pa^a^ete. . (a. we>. a. s.e A). He.ce, we do 
not show this case in Fig. IHl 

The variation of the binding energy with b in the three-body system and higher is a key 
indicator that V^i2,e// is only approximative for these systems. Furthermore, the fact that the 
calculated binding energies can overbind for certain values of b, as shown in the three-body 
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system, shows that this approximation is not variationaL For the alpha particle, our results 
underbind by ~ 1 MeV. 

B. P-shell nuclei 

Figure |3] shows our results for the ground state binding energies for ^He, ^Li, and ^Li 
usmg the Av8' potential 0- 

Again, our results are displayed with points, while the 'exact' 
theoretical binding energies are taken from Refs.Q] and displayed by solid lines. For the 
seven-body system, the calculated GFMC binding energy is -33.56 MeV, which is below the 
plotted range in Fig. |3J Numbers in parenthesis refer to the quantum numbers (J'^,T). For 
the five- and six-body systems, we show calculations for both A = and A = 2 included 
spaces. For the p-shell nuclei presented here our results underbind for all values of b. We plot 
our results where our calculated binding energies approach closest to the exact theoretical 
results. The selected ranges of b in Fig. |3] reflect this fact. Figure El shows the same results 
using the Avl8 potential. 

Since the A = calculations of these p-shell nuclei all underbind, they are susceptible 
to cluster breakup. Furthermore, their calculated binding energies are nearly degenerate. 
We do not have a satisfactory explanation for this behavior. This apparent degeneracy is 
broken in the A = 2 calculations of the five- and six-body systems. We have not performed 
any A = 2 calculations for the seven-body system due to time and computer limitations. 

P-shell nuclei have richer structure than their s-shell counterparts due to the existence 
of excited states. In Figs. |Bl|7|we show the excited spectrum of ^Li and ^Li using the Av8' 
potential. We also show the experimental spectrum to provide reference points, but stress 
that our results should only be compared to GFMC/NCSM results only. All BH calculations 
were done in A = included spaces. The quoted b values are taken at the minimum binding 
energy for each level. The AvlB results are shown in Figs. IHHni 

Our A = calculations of excited spectra fall within the range of GFMC and NCSM 
results. We note that the NCSM ^Li calculations were performed in a A=14 included space 
using HQ = 11 MeV (6=1.94 fm) and gave 28.08 MeV binding for the ground state of ^Li 
using Av8'. The excitation energies of the (2^,0) and, in particular, of the (1^,0) states 
were not converged yet and are still decreasing with increasing A. The NCSM ''Li results are 
taken from Ref. Q. In general our level orderings correlate with those of GFMC and NCSM 
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calculations as well. However, certain levels are inverted, such as the (3"'",0)-(0+,l) states 
of ^Li and the (3/2~,l/2)-(l/2^,l/2) states of "^Li. These particular levels are sensitive to 
the strength of the spin-orbit interaction in the A^A^ interaction, which to date, is poorly 
represented in all potential models. Clearly our approximation of Vejf in these limiting 
included spaces exacerbates this issue. The energy differences between using Avl8 and Av8' 
are typically less than ~ 200 KeV. 



VI. CONVERGENCE ISSUES 

The A = 2 calculated binding energies for the five- and six-body systems of Figs. EE 
are closer to the exact answers compared to the A = results, and the dependence on b 
is somewhat flatter. This suggests some sort of convergence with increasing A. Indeed, as 
A — > oo, one has that P — ^> 1, giving Vi2,e// V^. Hence in this limit the calculated binding 
energies are guaranteed to converge to the exact answers. However, we caution that since 
we are working in such small included spaces, such limiting behavior might not be seen until 
A becomes much larger. Furthermore, the argument above does not give any information on 
how quickly the binding energies converge to the exact answer. For example, larger model 
space calculations of the ^H fFig.lTUI) do not show a systematic convergence. This is in direct 
contrast to NCSM calculations, where there is definite convergence as the included space is 
increased jS]. A possible explanation for this difference is in the incorporation of a confining 
HO mean-field potential in NCSM calculations. The starting bare Hamiltonian for NCSM 
calculations is 

o-NCSM Tj , m I rrHO 
-n = U + Ic.m. + U^^, 

where H is given by Eq. ((2|) and Tcm. + U^m is the cm. harmonic oscillator Hamiltonian. 
In principle, the addition of a confining HO Hamiltonian (which is in the end subtracted 
from the effective Hamiltonian in any case) does not alter the final answer of the many-body 
problem, as long as the effective interaction is solved to a high level of accuracy. Indeed, 
we argue that the inclusion of a confining mean-field increases the rate of convergence as 
the size of the included-space is increased. Since our BH calculations lack any confining 
mean-field {e.g. U^^), our convergence with increasing A is hindered. 
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VII. CONCLUSION 



The exponentially increasing complexity of the many-body Hilbert space motivates the 
use of effective interactions in truncated included spaces. Unfortunately, the exact effective 
interaction -ffe// is as difficult to deal with as the original infinite-Hilbert space problem. 
Hence some form of approximation must be made to -ffe//, or more specifically, to Veff. 
In Sect, mil we showed a particular approximation of Veff that was obtained by summing 
ladder diagrams to all orders, giving Vi2,e//- We argued that our approximation of Veff, 
one in which spectator dependence and cluster recoil are manifestly included, was a natural 
extension of the single- valued G-matrix of traditional shell model calculations and the multi- 
valued G-matrix of Ref.[n|- Our approximation of Veff is equivalent to the leading order 
term of the BH expansion of Ref . [lO| ■ 

Because we include the effects of cluster recoil, our approximation of is an A-body 
operator, as opposed to the two-body nature of the G-matrix interactions mentioned above. 
In Sect. IIVI and Appendix El we showed how A-body matrix elements of Vi2,e// can be 
easily calculated by transforming the A-dimensional nested integrals into a sum over one- 
dimensional integrals. The result is valid for any A-body system and can be computed 
quickly and accurately. We showed that by integrating over the spectator degrees of freedom, 
we could define an average residual interaction for the interacting two-nucleon cluster that 
depended on the initial HO momentum configurations of the spectator nucleons. Figure |21 
shows how this average interaction in the ^Sq channel changes as the number of spectator 
nucleons is increased. 

Finally, in Sect.|V|we showed our results of using V\2,eff in s- and some p-shell nuclei for 
the smallest allowed included spaces. All our calculations produced stationary solutions [i.e. 
negative binding energies). Our s-shell results show good agreement with exact theoretical 
answers. In particular, the alpha particle is underbound by only ~ 1 MeV. On the other 
hand, our p-shell results all severely underbind in A = included spaces, and are therefore 
susceptible to cluster breakup. The excitation spectra of ^Li and ^Li fell within the ranges 
calculated by GFMC and NCSM methods. The level orderings correlated to some extent 
with these methods as well, though level inversions do exist. We argued that these level 
inversions were closely related to the weak spin-orbit interactions of the potential models 
used in our calculations, which were further exacerbated by our approximation of Veff in 
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such small included spaces. 

For both ^He and ^Li, we extended our calculations to A = 2 included spaces (see Figs.Ol- 
E}, and some semblance of convergence was observed. We caution, however, that any appar- 
ent convergence in such small spaces is not definitive. Indeed, for the case of '^H, convergence 
was not seen as the size of the included-space increased, as shown in Fig. [Ul This is in direct 
contrast to NCSM calculations. Though our calculations are, in a sense, more sophisticated 
than NCSM calculations, our results are no better for A > 2. We argued in Sec. IVII that the 
lack of a confining 'mean-field' potential was the most likely reason for our poor convergence. 

As mentioned in the previous section, our calculations were not too computationally de- 
manding since we worked within the smallest allowed included spaces. Indeed, we were able 
to perform all our calculations on a single Unix box. Calculations on higher A-body systems, 
or on much larger included spaces, will require that our calculations become parallelized, 
which we hope to do in the near future. It would be interesting to see if the A = 'de- 
generacy' persists for larger p-shell nuclei. To tackle our convergence issues, we intend to 
incorporate a 'mean-field' potential into our future calculations. 



APPENDIX A: A-BODY MATRIX ELEMENTS OF Z 



eff 



The relevant expression in Jacobi coordinates is 



X 



(Al) 



where Rni{p) is the HO radial wavefunction in momentum space. Replacing the radial 
integrals by their series expansion and changing to dimensionless variables xi, . . . ,xa-i, 
Eq. lAll becomes 



mi,m' ■•• ,mA-i,m'._,=0 



X 



dxi ■ ■ ■ dxA-ie' 



2{h+mi+m[+l) 2(«A-i+mA-i+"i^_i+l) 



■ X 



A-1 



E 



(A2) 
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where 



'A-1 



n 



^2T{n, + l)r(n, + h + 3/2)) 
\ r{ni-mi + l)r{h + mi + 3/2)r{mi + l) 



'A-l 

n 



v/2rK + i)rK + /, + 3/2)) 



• (A3) 



A r«-m^ + i)r(/i + m^ + 3/2)rK + i) 

The integral in Eq. lAll can be evaluated by first changing variables to (A-l)-dimensional 
polar coordinates, p, 9a-5, ■ ■ ■ ,0i,4>, where < < 27r, and < 9i < n: 

Xi = p coscj) sinOi sin62 ■ ■ ■ sin9A-5 
X2 = p sincj) sinOi sin92 ■ ■ ■ sin6A-5 
xs = p cos 9i sin92 ■ ■ ■ sin9A-5 



Xi = p cos9i_2 sin9i_i ■ ■ ■ sin9A-5 



(A4) 



xa-1 = p cos9 



P 



A-5 



E 



The infinitesimal volume element is 



dxi ■ ■ ■ dxA-i = p"^ sin^ ^9a-5 ■ ■ ■ sin9i dp d9A~5 ■ ■ ■ d9i 



With these change of variables, the integral in Eg. lAll becomes 



tt/2 



X 



tt/2 



tt/2 



(cOSeA-5)'^'^-^+'"^-^'""^-^+'Hs^n^A-5)'^^="^'^'^+'"^+'"^+'^'"^~'rf^A-5 



X 



dp e 



, pA-2+2 Y,fS^\k+m,+mr+l) 



E- 



• (A5) 
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Note the limits of integration. All the angular integrals and the integral over p can be 
expressed analjd;ically, 

f7r/2 



dd cos'^9 sin^ 



2a+l 

dp e-f^ 



r(^)r(^) 



^ 2 



2r(^^±^ + i) 



'0 e-p^ 2 

Plugging Eqs. IA5I - 1 A 71 into Eq. IA2I gives the final result 



e"e-' r(-a,-e)r(a + l). 



(A6) 
(A7) 



/ 7/ / 7/ 



1 



E-T 



{nih ■ ■ -ua-iIa-i >- 



ni,n, ,ra^_i,n , 



E 





1=0 




r2E\ 







huj' 



(A8) 



where a = (A — 3)/2 + + mj + + 1) and 



E(mi, m^, ■ ■ ■ , mA-i, m^_i) 



\Si{ni+mi) 



'A-l 



n 



v/2r(n, + l)r(n, + + 3/2)) 

\ r(ni-m, + i)r(/i + mi + 3/2)rK + i) 

V2r(n^ + l)rK + /, + 3/2)) 



n 



\ r«-m^ + i)r(/i + m^ + 3/2)rK + i) 



n 

i=2 



T{li + mi + m'i + 3/2)r(^* + m^- + m'^ + 1) + 
2r(/, + mi + + Etl(^, + + m' + 1) + ^) 



■ (A9) 



APPENDIX B: INTEGRATING OUT SPECTATOR DEGREES OF FREEDOM 



The relevant expression here is 

< p'] n[l[, ■ ■ ■ , ^A-2^A-2l-^b; ^ih, ■ ■ ■ , nA-2lA-2 >, 



(Bl) 



where the interacting two-nucleon initial and final momenta are p and p', respectively, and 
the spectator quantum numbers are denoted with subscripts that run from 1 through A — 2.. 
F is some operator that depends only on p, p', and where qi are the spectator 

momenta. Gq, ti2, and Go^i2G'o are three examples of such an operator. Since F only 
depends on qf, which is a scalar, it cannot change the orbital quantum numbers of the 
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spectator nucleons. Hence Eg. IB II becomes 



X J dqi--- dqA-2qf ■ ■ ■ ql-2K{h (^O^m^i (^i) " " " K'^_^Ia-2 {<lA-2)RnA-2iA-2 (^^-2) 

A-2 

X F{p,p'-Y,qf)- (B2) 



i=l 



Equation IB2I is very similar to Eg. lAll Indeed, we can follow a similar procedure outlined in 
the previous section to come to our desired result. Repeating the steps from Eg. I A2I through 
Eg. IA7I but with a change of variables to (A-2)-dimensional polar coordinates rather than 
(A-1), gives 

Si[,h ■ ■ ■ ^i'a-2^1a-2 Y1 C{mi, m[,--- , mA-2, m'^-2) 

Jo " J 

where b is the oscillator parameter and 



C{mi,m[, ■ ■ ■ , mA-2, m'j^-2) 



'A-2 



n 



^/2T{n, + l)r(n, + + 3/2)) 
\ r{ni-mi + l)r{h + mi + 3/2)r{mi + l) 

yr v/2r(n- + l)rK + /. + 3/2)) 

11 r«-m^ + i)r(/, + m: + 3/2)rK + i) 



'A-2 

n 

.i=2 



T{k + mi + m[ + 3/2)T{J^-J^{lj + ruj + m'j + 1) + 



i+2^ 



(B4) 



Note that the integral in Eg. lB3l cannot be determined analytically, as in the previous section. 
However, it can be numerically evaluated guickly and accurately. 
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FIG. 1: (Color online) Diagrammatic terms kept in ladder approximation of Eq.EI Nucleons 1 and 
2 are interacting via a ladder sum in the presence of spectator nucleons 3-A. Dashed lines represent 
interactions with Vi2- 



TABLE I: Dimensions of A = 0, 2 included spaces for different few-body ground states. 



A 


2H 


^H/^He 


^He 


^He 


6Li 


^Li 





1 


1 


1 


1 


3 


5 


2 


3 


5 


5 


26 


48 


185 



20 




FIG. 2: (Color online) Spectator dependence of Goti2G'o. All axes are dimensionless. Plot (a) is 
calculated with no spectator nucleons. Plot (b) includes two extra spectator nucleons on top of the 
interacting pair. In this case, the Jacobi relative angular momentum for these nucleons are in the 
Os states. Plot (c) and (d) include a third and fourth nucleon, respectively, in the Jacobi relative 
Op states. 
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FIG. 3: (Color online) S-shell results for A = included spaces. Points cor resp ond to our calcula- 
tions, whereas solid lines correspond to the 'exact' results taken from Ref. la]- In the three-body 
system, the upper and lower set of points refer to ^He and '^H, respectively. We do not show results 
for ^H, since all calculated points lie on the exact theoretical result, regardless of b. Calculations 
were done with the Avl8 potential. 



22 



He 



> 




-5 



-10 



-15 



-20 



-23 . 85 
-25 



-30 



-1 1 1 1 r 



(3/2 ,1/2) 



'Li 



A=0 



A=2 

xxxxxxxx 



_] I I I L_ 



> 
d) 



-5 



-10 



-15 



-20 



-25 



1.3 1.35 1.4 1.45 1.5 1.55 1.6 
b (fm) 



-28.19 
-30 



-1 1 1 1 r 



(1,0) 



A=0 



A=2 



X X X X 



_] I I I L_ 



> 




-10 



-15 



-20 



-25 



1.5 1.55 1.6 1.65 1.7 1.75 1. 
b (fm) 



-30 



-1 1 1 1 r 



(3/2 ,1/2) 



A=0 



_l I I I L_ 



1.6 1.65 1.7 1.75 1.8 1.85 1.9 
b (fm) 



FIG. 4: (Color online) Same as in Fig. Bl but for p-shell ground states using Av8'. Solid lines 



correspond to results taken from Ref. 12 
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FIG. 5: (Color online) Same as in Fig. \^ but using Avl8. 
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FIG. 6: (Color online) Excited spectrum of ^Li using Av8'. 
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FIG. 7: (Color online) Excited spectrum of ^Li using Av8'. 
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FIG. 8: (Color online) Excited spectrum of ®Li using AvlS. 
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FIG. 9: (Color online) Excited spectrum of ^Li using Avl8. 
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FIG. 10: (Color online) Dependence of '^H on models space size A and oscillator parameter b using 
Avl8. The solid line is the exact '^H binding energy with the Avl8 potential. 
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